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We investigate non-Gaussian states of light as ancillary inputs for generating nonlinear transfor- 
mations required for universal quantum computing with continuous variables. We consider a recent 
proposal for preparing a cubic phase state, find the exact form of the prepared state and perform 
a detailed comparison to the ideal cubic phase state. We thereby identify the main challenges to 
preparing an ideal cubic phase state and describe the gates implemented with the non-ideal pre- 
pared state. We also find the general set of gates that can be implemented with ancilla Fock states, 
together with Gaussian input states, linear optics and squeezing transformations, and homodyne 
detection with feed forward. Such circuits can be used to approximate a certain class of Hamil- 
tonian evolutions. Furthermore, we analyze the question of efficient classical simulation of these 
gates. These results extend the existing theorems about efficient classical simulation for continuous 
variable quantum information processing. 



Universal quantum computation can be performed by encoding qubits in discrete states of quantum systems and 
performing a set of universal quantum logic gates or unitary operations on the qubits. For the case of light, the 
electromagnetic field amplitudes are described by operators that have a continuous rather than discrete spectrum. 
Quantum information processing with continuous variables (CV) is thus appealing because of the possibility of optical 
implementation using simple linear optics elements and coherent states of light. Furthermore, a CV quantum computer 
may be more suitable for simulating CV quantum systems. Whereas much progress has been made in areas of 
communication, measurement and computation using optical systems jT], [H, universal CV quantum computation 
with light remains a challenge due to the weakness of optical nonlinearities. This problem can be circumvented by 
off-line preparation of nonclassical states of light that are fed into the optical circuit at various points to generate 
the required strong nonlinearities. This idea was explored by Knill, Laflamme and Milburn [3] and subsequently by 
others [1, [1] for quantum computing with qubits, and by Gottesman, Kitaev, and Preskill (GKP) for qudits encoded 
in CV systems Q ■ In this paper we present an analysis of non-Gaussian ancilla states of light in optical circuits as 
an off-line resource for generating universal gates for continuous variable quantum information processing. 

For universal CV quantum information processing, the aim is to apply a finite set of gates that, by repeated 
applications, will generate unitary operations corresponding to a general class of Hamiltonian evolution. Specifically, 
we are interested in Hamiltonians that are arbitrary polynomial functions of the continuous variables [sl]. In optical 
implementations, combinations of linear optics and squeezing elements can generate transformations that are limited 
to quadratic functions of the field amplitude operators x and p. Furthermore, any circuit consisting of input Gaussian 
states and these transformations, along with homodyne detection and feed forward can be efficiently simulated on a 
classical computer Q. This efficient simulation is due to the fact that the process can be described as a Gaussian- 
preserving map and can be completely characterized by the evolution of the means and covariances of the initial 
Gaussian state. This requires performing simple operations on a finite set of variables, that scales polynomially with 



Nonlinear transformations, namely those outside the set of linear and quadratic operations, are therefore required 
for universal CV quantum computing - i.e. for implementing unitaries corresponding to Hamiltonians that are poly- 
nomial functions of the continuous variables. One possibility is to use photon counting measurements to generate 
nonlinearities [TT| . Photon counting plays a key role in schemes for linear optics quantum computation using qubits [1| 
and for computation with CV-encoded qudits The measurements can be used to generate non-Gaussian states 
through an 'off-line' preparation step that is independent of the particular circuit being implemented. GKP then 
showed that these non-Gaussian states can be employed as ancillae at different stages of 'on-line' quantum informa- 
tion processing to replace the requirement for nonlinear gates. This partition of the quantum computation into off-line 
and on-line elements obviates the requirement for real-time photon-counting measurements that must adapt to the 
vagaries of each specific circuit. 

Here we investigate and analyze off-line preparation of non-Gaussian ancilla states and the employment of such 
ancilla states in on-line circuits. We specifically consider cubic phase states and Fock states. Our interest focuses 
on these off-line prepared states as resources for replacing the nonlinearity in order to effect universal quantum 
information processing. We first perform a detailed analysis of the scheme proposed by GKP for off-line preparation 
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of a cubic phase state. We find the exact form of the state generated by the GKP scheme and show that it is a finite 
superposition of displaced number states. We use the fidehty as a figure of merit to assess the efficacy of preparing 
an ideal cubic phase state with this scheme, and present detailed studies of the Wigner function for experimental 
parameters that take into account finite squeezing and finite detector efficiency. Our analysis unfortunately indicates 
that current levels of squeezing arc insufficient to prepare close approximations to the cubic phase state, although 
this limited squeezing does not necessarily rule out this approximate cubic phase state from providing an off-line 
resource for universal quantum information processing with Gaussian operations in the on-line circuit. Our results 
are of importance to current experimental research with squeezed light and in particular, it provides guidelines to 
those actively working on experimentally preparing the cubic phase sta te Il2| . This work is also of relevance to recent 
studies of quantum computing with continuous variable cluster states [l3j . 

We also consider an alternative non-Gaussian state as an off-line resource for universal quantum information pro- 
cessing with Gaussian elements in the on-line circuit, namely the ubiquitous Fock state. Whereas the cubic phase 
state is special in that it both delivers an off-line resource for universal quantum information processing with Gaussian 
circuit elements and also makes the incorporation of the cubic phase state into the circuit simple and elegant, the Fock 
state is also non-Gaussian and perhaps the most feasible to prepare, at least for single-photon states [14]. Since the 
Fock state is probably the most amenable non-Gaussian state to feed into a Gaussian circuit as an off-line resource, 
it is important to assess its status as a useful off-line resource for CV computation. 

In this analysis of Fock state resources, we consider a Gaussian-preserving circuit receiving ancillary photons. We 
determine the general form of non-Gaussian output states that can be obtained from a Gaussian-preserving circuit 
with an n-channel input state along with some ancillary channels, each occupied by a number state. Our studies show 
that ancilla Fock states appear to provide the requisite nonlinearity for universal CV computation, but circuit design 
for implementing nonlinear gates is not obvious with Fock state inputs. We also include a brief discussion of efficient 
classical simulation of such maps by considering the classical resources required to simulate these maps. Our results 
suggest that such maps may not be possible to simulate efficiently classically, as one would expect for non-Gaussian 
circuits. 



II. THE CUBIC PHASE STATE ANCILLA 



A. The cubic phase gate 

Gottesman, Kitaev, and Preskilll (GKP) define the (unnormalized) cubic phase state as Q 

\^)= fdxe'^^'lx) (1) 



and propose its employment as an off-line resource for the Gaussian-preserving circuit depicted in Fig. 1. Here the 
input state is in a position (x) eigenstate, which is a limiting case of displaced squeezed states, with the squeezed 
state given by 

h)=exp(^|at2-!^a2)|0). (2) 

The state \x) is approximated for 77 G M becoming very large and displaced by the unitary displacement operator 

D{a) ~ exp(aa^ — a*a) (3) 

with a = a; e M. The cubic phase state, I7) serves as the ancilla. The Gaussian-preserving sum-inverse (SUM~^) 
gate is defined by the inverse of SUM = exp(— ii;ip2)j where the subscripts label the input and the ancilla channels 
respectively. 

The importance of the cubic phase state is that it can serve as an off-line resource for effecting a cubic phase gate, 
namely exp(i7x'^), on the input x eigenstate using only the cubic phase state, the Gaussian-preserving sum-inverse 
gate, a homodyne measurement of the ancilla and feedforward. After applying the SUM~^ gate and measuring the 
ancilla channel in the x basis to obtain outcome q, the resultant output state is 



IV'out) = (9I SUM-1 la:) I7) 

= J da;2e'^"^ (q| SUM^^ |x) {x^) 



e '^ X) 
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FIG. 1; Circuit for implementing a cubic phase gate on an input x eigenstate using an ancilla cubic phase state, and Gaussian- 
preserving gates. The gate U{q) depends on a measurement of the ancilla in the i-basis with outcome q. 

Then the cubic phase gate can be effected by applying a quadratic (i.e. Gaussian-preserving) g-dependent Hamiltonian 

H = iq-l{x'^ + qx + q^/i), (5) 

which generates the unitary transformation 



?7(q)==exp|i — 7 (x -I- gr)^ | 



(6) 



This unitary operator acting on the output state (4) yields e''''^ |a;) 



B. GKP scheme for preparation of the cubic phase state 

The cubic phase gate, plus the application of linear optics and squeezing operations, homodyne measurements and 
feed forward, is sufficient to generate a set of gates for universal CV quantum information processing by optical means. 
However the ideal cubic phase state defined in Eq. (1) is not normalizable and therefore does not correspond to a 
physical state. GKP address this issue by proposing a scheme, depicted in Fig. 2, to prepare an approximation to 
the cubic phase state via displacement and photon counting on one channel, or mode, of a two-mode squeezed state; 
the other mode approximates a cubic phase state with a parameter conditioned by the photon counting measure- 
ment outcome on the other state. The two-mode squeezed state is prepared from the vacuum state by the unitary 
transformation 



S,j{z) 



exp<; ^a\a] ~ ^a^a^ 



(7) 



analogous to the single-mode squeezed state explained in Eq. We analyze the scheme in detail here. 

Consider the circuit shown in Fig. 2. We start with two input vacuum states, which then undergo two-mode 
squeezing 6*12 characterized by a real squeezing parameter z = 2r S M in Eq. ([7]). One of the two states is mixed 
with a coherent light field at a beam splitter, which results in a momentum displacement D2{a2 = ipo) for D given 
by Eq. ©. 

This channel is subject to photon counting with measurement outcome n2. Initially we consider the ideal case of 
perfectly efficient photon counting. For the case of ideal measurement, the state in the other channel conditioned on 
this measurement can be written in the x representation as 



= {xi,n2\D2{a2)Si2{r) |0,0) 



(8) 



For the input state expressed in the x representation, we have 



%l]'('"^\xi) ^ I da;2 exp(ia2X2)/(a;i,a:2)M„2(a;2), 



(9) 

with M„2(x2) the Hermite Gaussian functions [l5| . Here f{xi,X2) is the x- representation of the two- mode squeezed 
state 



fixi,X2) =iVexp{ - ^[coah{2r){xl + xl) 
+ 2 sinh(2r)xia;2]} 

with N a normalization constant. 

GKP approximate the integral in Eq. (9) in the large r, large 712 limit as 



(10) 
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FIG. 2: Circuit for preparing a cubic phase state 
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C. Exact analysis of GKP preparation scheme 



(12) 



We compute an exact expression for the state 'ip''^^\xi) prepared via the GKP scheme (for convenience we caU this 
the GKP state). We use the generating function for the Hermite Gaussian (see Appendix) so that 



V'^^(xi) 



gn 



dx2 e2^^*-^'/V"^^V(a;i,a;2) 



(13) 



As f{xi,X2) is a Gaussian function of X2-, the above integral can now be calculated. After some manipulation, the 
rij'^-order derivative of this integral results in an expression for the state given by 



with 



gniix) 



^ ) (— tanhr)™ (iposech^r)"^ ™Um(a;). 



(14) 



(15) 



D. Comparison between GKP state and ideal cubic phase state 

To scrutinize the complicated state in Eq. (14), we first examine the m = n2 term of the sum in Eq. (14). For large 
712, far from the turning points, we can use a semi-classical approximation for Un2{x) (see Appendix). Then, close to 
X = 0, if posech^r <C 1 the m — n2 term dominates and the state can then be approximately written as 

i/'("^)(xi) oce'^^^? +e'^^i (16) 



with 7 given by Eq. (12) and (i — 2^2n2 + 1. 

The required cubic phase term thus arises from the m = n2 term in Eq. (14). When other terms in the sum are 
included, they can act to reduce the contributions of the extra linear phase term in the large r large n2 limit, thus 
leading to an approximate cubic phase state. We perform a numerical analysis to study the conditions for which the 
prepared state approaches a cubic phase state by computing the corresponding Wigner function [T?! given by 



W{x,p)= J dyip* {x — y)ipix + y)e' 
The Wigner function of the ideal cubic phase state is 



-2ipy 





4 


1/3 / 


■ 4 " 


2ttN 




A. 






3^ 




.37. 



1/3 



[3-fx^ 



(17) 



(18) 



5 



(a) W 




-1 1 2 3 4 5 
P 



FIG. 3: (a) Wigner function of the ideal cubic phase state for 7 = 0.05. (b) Overhead view, (c) a; = shce. 



with Ai(a;) the Airy function and N a normahzation factor. The Airy function causes the Wigner function to oscillate, 
resulting in 0-dependent negative regions with e = tan-i(p/a;) (Fig. 3). 

The Wigner function of the prepared state described by Eq. (14) in the regime of large squeezing and large 712 is 
a sum of the Airy function of Eq. (18) obtained from the m = n2 term and additional functions arising from the 
additional terms in Eq. (14). In order to examine the contribution of these terms for a finite amount of squeezing, 
we compute the Wigner function of the prepared state. Fig. 4 shows the effect of the squeezing parameter on the 
negativity of the Wigner function of the prepared state along the x = line for 71,2 w 50. The data show that, 
for currently achievable levels of entanglement [H, [l^, HO, [HI , the oscillations in the Wigner function are strongly 
suppressed and the state does not closely approximate a cubic phase state. Furthermore, the analysis thus far assumed 
perfect efficiency of the photodetector in channel 2 in the circuit of Fig. 2. In reality one must take into account an 
imperfect photodetector. 

If the photodetector has an efficiency 77, then the output state conditioned on a measurement of 712 photons in 
channel 2 is a mixed state that can be written as [22j 

00 

p^^J2 Pn\4"^){4"\ (19) 
n=n2 

with 

Pn - • ,, V''' (1 - ?7)("-"^) (20) 

n\[n — 7^2)! 

the BernoulUi distribution. Fig. 5 shows that as the detector efficiency decreases, the oscillations in the Wigner 
function are further washed out so that the actual state diverges from the ideal cubic phase state. 

The threshold theorem of quantum error correction requires that quantum gates be implemented with a minimum 
fidelity ^23]. If we use the prepared state given by Eq. (16) as the ancilla in the circuit of Fig. 1, we obtain a 
superposition of a cubic phase gate and a quadratic transformation on the input state. In general, the non-ideal 
prepared ancilla state described by Eq. (14) will result in a gate that is a superposition of the desired cubic phase 
gate and additional gates. The fidelity of implementing a cubic phase gate will be limited by the fidelity of preparing 
the cubic phase state I7). The ideal cubic phase state is not normalizable, making it challenging to assess the fidelity 
of preparing this state. Here we define the fidelity normalized over the domain D ranging from Xmin to Xmax, where 
the prepared state p{x,x') = {x\p\x') has non-negligible support. For numerical evaluation, we define 

^'^^'^''"1 otherwise " ^ > 
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FIG. 4: Wigner function of the prepared state for different values of the squeezing parameter. 



The fidelity F is then calculated to be 



with 



F : 



D = 



D ' 



dx\{x\-f)\ 



(22) 



(23) 



We find that for a squeezing parameter r = 1.34 and perfect detector efficiency, the normalized fidelity is roughly 
20%. Hence, our studies show that a close approximation to the cubic phase state cannot be prepared using current 
levels of squeezing and detector efficiency. 

We can define a corresponding fidelity of performing the cubic phase gate using the circuit in Fig. 2 for a given 
input state lip) as 



F, 



gate 



= Y^(^|[/taC/|^). 



Here a = U\ip){^\U\ 



(24) 



(25) 



and 



UH) = U{q) {q\ SUM-l IV) 



(26) 

with U{q) given by Eq. (6) and being the prepared state given by Eq. (14). The minimum gate fidelity is the gate 
fidelity minimized over all input states. Writing the input state in the x basis, the gate fidelity becomes 



- gate 



dx |?A(a;)pe-'^("+«'V^"^^(a; + g) 



(27) 



with = and ^p[^^\x) given by Eq. (14). For an input state that is an equal superposition of x eigenstates, 

the state fidelity in Eq. (22) is equivalent to the gate fidelity for q values close to 0. Hence, as one would expect, the 
fidelity of preparation of the ancilla state constrains the minimum gate fidelity of the cubic phase operation. 

We note that the GKP scheme for preparing the cubic phase state uses the minimal number of non-linear elements 
- the key element which generates a non-linear operation is a single photodetector. Any other scheme that relied on 
Gaussian operations and required more photodetectors with finite efficiency may reduce the performance. However, 
the use of non-Gaussian ancillae in the preparation scheme may improve the preparation fidelity as discussed in the 
following section. 

Whether or not the non-ideal prepared state can be used efficiently for universal quantum information processing 
with Gaussian preserving maps is an open question; perhaps a scheme using the prepared state in the circuit in Fig. 
1 together with some form of error correction could provide a solution to this problem. A possible alternative might 
be to prepare many copies of the non-ideal state and perform some purification protocol in order to prepare a state 
closer to the ideal cubic phase state. Furthermore, the use of alternative schemes using other non-Gaussian ancilla 
states such as Fock states may generate better approximations to the cubic phase state. We discuss this further in 
the following section on Fock states. 
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FIG. 5: Effect of finite detector efficiency on the Wigner function of the prepared state for r = 1.34. 

III. THE FOCK STATE ANCILLA 
A. On-line operations generated with ancilla Fock states 

Given the difficulty of preparing a cubic phase state via the GKP scheme, we explore the use of other non-Gaussian 
states as possible ancilla states in Gaussian-preserving circuits for implementing universal CV computation. To this 
end, we analyze the transformations generated by input ancilla number (or Fock) states. We consider a general optical 
circuit as shown in Fig. [6Ka). The input state |i/'in) consists of n electromagnetic field modes, each described as a 
harmonic oscillator. The field amplitude operators x and p for each mode obey the commutation relation [x,p] = ih. 
We assume that the input state is a product of Gaussian states for each of the n channels. Eigenstates of x or p can 
be regarded as infinitely squeezed Gaussian states and are included in the set of possible input states. 

We first analyze the gates generated by a single ancilla Fock state |n2) with a specific number of photons n2. The 
operator U12 is any element of the Clifford group for continuous variables Q constructed from a combination of linear 
optics elements and squeezers. The output state is conditioned by a homodyne measurement of the ancilla, such that 

IV'out) = (y2|f/l2|^in) 1^2). (28) 

Here subscript 1 labels the input state, and subscript 2 labels the ancilla. Rewriting the number state \n2) as a 
derivative of a coherent state \a2{t)) with 02 (^) = V^t (see Appendix), 

(29) 

Inserting this into Eq. (28), we obtain the output state 

(30) 
where 




\n2) = j^exp{f)\a2{t)) 



Qn2 

|V'out) = 7^exp(t2)|x(0) 



lx(i)) = (y2|c^i2lV'i„)l«2(i)) 



(31) 



The state |x(i)) is a Gaussian state since it results from a Gaussian-preserving map [l0'| with input Gaussian states 
undergoing CV Clifford group transformations followed by homodyne measurements. The output state in Eq. (30) can 
also be derived by using the result that any linear optics circuit can be decomposed into an interferometer, followed 
by a parallel set of single mode squeezers, followed by another interferometer and additional displacements [24j . 

We now consider adding an ancilla prepared in a product of number states rather than a single number state. 
Fig. 6(b)) presents the case of N photons entering the circuit in fc — 1 channels with Ya=2 = Then, using the 
expression for a number state given by Eq. (29), we obtain a general expression for the output state at the end of the 
computation, 
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FIG. 6: (a) Gaussian-preserving circuit for using a single Fock state with n2 photons as an ancilla. A homodyne measurement is 
performed on the ancilla with outcome j/2 • (b) Gaussian-preserving circuit for using A'' photons distributed over k — 1 channels 
as ancillae. The circuit is conditioned on homodyne measurements performed on the ancillae with outcomes 1/2 ■■ - yk- 



where 

\x{t2,-tk)) = Uik{yk\Uik...{y3\Ui3U{y2){y2\Ui2\^in)<E)\a2{t2)) ...\(^k{tk)) (33) 

is a Gaussian state resulting from Gaussian-preserving unitaries and homodyne measurements acting on the input 
Gaussian state. In the a:-basis, the derivative operations in Eq. (32) result in an output state which is a product of a 
Gaussian function and a polynomial function of x of maximum order 2(fc — 1) for fc — 1 ancilla channels. By changing 
the circuit parameters, different such functions can be generated. 



B. CV computation with ancilla Fock states 



It is not immediately obvious from the form of the output state in Eq. (32), whether this circuit leads to universal 
CV quantum computation or not. In order to analyze whether this circuit is sufficient for universal CV quantum 
computatioii, we only need to consider a single-mode rather than a multi-mode input. As shown by Lloyd and 
Braunstein [Si, any non-Gaussian unitary operation on a single mode together with the set of Gaussian preserving 
operations is sufficient for universal CV quantum computation. This implies that since the output states in Eq (32) 
are not Gaussians states, the circuit in Fig. 6 can generate a non-Gaussian operation on a single input mode and 
hence appears to be sufficient for CV quantum computation. The challenge then is to design a particular circuit 
that can generate a useful nonlinear (non-Gaussian) operation, for example a cubic phase gate on a single mode. In 
fact, our analysis of the GKP preparation scheme for the cubic phase state allows us to identify one possible simple 
circuit to perform a non-Gaussian operation on a single mode using ancilla Fock states. Consider a single mode input 
state and a second ancilla mode having a single ancilla Fock state with photon number 712. The ancilla state is first 
acted on by a momentum displacement operator with the displacement p — \/2n2 + 1. For large 712, the resulting 
displaced number state can be written in the x-basis close to a; = as a superposition of a cubic phase state and a 
linear phase state as given in Eq. (16). Thus such an ancilla state when fed into the circuit of Fig. 1 will generate a 
superposition of a cubic phase operation and a Gaussian operation on a single mode. Although this is not an ideal 
cubic phase operation, as noted before, together with error correction or purification, it may be sufficient for useful 
CV computation. 

We identified the circuit discussed above by noting that the displaced ancilla number state is none other than the 
TO = 712 term in Eq. (14) and that the state given in Eq. (14) is a finite superposition of Fock states that is displaced 
in momentum. Thus Eq. (14) shows that a finite superposition of displaced Fock states rather than a single displaced 
Fock state will closely approximate a cubic phase state for large 712 and r —^ 00, and hence generate a cubic phase gate 
when fed into the circuit of Fig. 1. Although the GKP scheme for preparing this state is limited by finite squeezing, 
other schemes for preparing arbitrary finite superpositions of Fock states [2^, [2^, [13, US [H, Ho'l may be able to 
overcome the squeezing constraint. Whereas optical schemes for preparing an arbitrary, finite superposition of Fock 
states are non-deterministic [25l . [26j , the interaction of light with atoms in a cavity can deterniinistically generate 
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arbitrary single mode optical fields [13, [H, [3 ^'^^ even Fock states with large photon numbers [s^ . The feasibility 
of preparing a good approximation to the cubic phase state via such schemes is the topic of future work. 



C. Analysis of classical simulation of Gaussian circuits with ancilla Fock states 



It is widely believed that universal quantum computation cannot be classically simulated efRciently. In order to 
further test whether the circuits in Fig. 6 can perform universal CV quantum computation, we perform a natural 
adaptation of the efhcient classical algorithm for simulating evolution when the circuit inputs (including ancillae) are 
all Gaussian states [9] to the case of Fock state ancillae. We find that the inclusion of Fock state ancillae renders 
this algorithm inefficient, which would be consistent with the conditions for universal quantum computation. As 
shown in Fig. 6 we look at the case of n input channels initially in a product of Gaussian states along with N total 
photons shared in a product state among k—1 ancilla channels. From [9|, we know that the Gaussian-preserving map 
that transforms the input state \ipin) into the state |x) in Eq. (28) can be efficiently simulated classically by simply 
following the evolution of the means and covariance matrix of the input state. For n input channels this involves 
tracking a polynomial number O(n^) of real values. In the computational (x) basis, the state |x) is thus a Gaussian 
function whose coefficients are easily computed from the means and covariances. 

Assuming that the subsequent derivative operations can be performed efficiently, we seek to find the resources 
required to characterize the final output state after all the derivative operations have been performed. The derivative 
operations result in a product of the Gaussian function and a polynomial function of order 2N for N ancilla photons. 
It is simple to see that the maximum number of coefficients required to specify a general polynomial of order 27V in 
n variables (for n input channels), is given by 



Numerical calculations verify that this quantity scales exponentially with the number of ancilla photons and input 
modes. 

Our results indicate that Fock states may be a useful resource for CV quantum computation. The goal of efficiently 
creating Fock states on demand even at the level of single photons is thus of primary importance both for discrete as 
well as CV quantum information processing with linear optics. However, the design of simple gates with Fock state 
ancillae is not obvious so even if Fock state ancillae provided a resource for universal quantum computation, circuit 
design remains a challenge. 



We have analyzed the off-line preparation and use of non-Gaussian states of light as ancillary resources in optical 
circuits for generating nonlinear gates for universal quantum computation. Our studies of the scheme to practically 
prepare a cubic phase state reveal that a good approximation to the ideal cubic phase state is not achievable with the 
best known current levels of squeezing and photodetectors. Nevertheless, the non-ideal prepared state can generate 
interesting superpositions of Gaussian and nonlinear gates that, together with error correction protocols and purifi- 
cation techniques, may be sufficient for universal quantum information processing. The question of designing more 
practical schemes for preparing and using this state remains open and requires further investigation. 

We have also performed a detailed analysis of ancilla Fock states in optical circuits, and obtained the general 
set of gates that can be implemented with such states. Our results indicate that a circuit using input Gaussian 
states, ancilla Fock states and consisting of Clifford group elements, homodyne measurements and feed forward could 
generate the nonlinearity required for CV computation, but design of circuits for performing nonlinear gates is not 
trivial, and needs further investigation. We also find that a classical simulation of the states generated with an ancilla 
which is a product of ancilla Fock states can scale poorly with the number of ancilla photons and input modes. 
This result extends previous CV quantum information no-go theorems 9, 10] and implies that further analysis of the 
transformations accessible with ancilla Fock states is of interest to assess their use for specific information processing 
protocols. 




(34) 



IV. CONCLUSIONS 
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APPENDIX 



APPENDIX: NUMBER STATES IN TERMS OF GAUSSIAN STATES 



The number state \n) can be written in the x-basis as 



\n) = / dxun{x) \x) , 



(A.l) 

where m„(x) is a Hermite Gaussian function describing the n^^ eigenfunction of a harmonic oscillator [15]. Using the 
relation 



Unix) = e ' T^e 



+2xt 



(A.2) 



t=o 



and interchanging the order of the integral and the derivative, the above equation can be rewritten as 

\n) = ^e-^' j Axe'^'--'/'\x) 



exp(r) / da;e 2 



(A.3) 



Recognizing that the integral is a representation of a coherent state \a{t)) with a{t) = ^/2t, we can express the number 
state as 



l'^) = ^exp(i2)|a(t)) 



For large n, far from the turning points, we can use a semi-classical approximation for Un{x) to write 

^2T:pn{x) 



(A.4) 



(A.5) 



with Pn{x) — \^2n + 1 — x"^ 
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